Disproportionation and Metallization at Low-Spin to High-Spin Transition in 

Multiorbital Mott Systems 
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We study the thermally driven spin state transition in a two-orbital Hubbard model with crystal 
field splitting, which provides a minimal description of the physics of LaCoOa. We employ the 
dynamical mean-field theory with quantum Monte-Carlo impurity solver. At intermediate temper- 
atures we find a spin disproportionated phase characterized by checkerboard order of sites with 
small and large spin moments. The high temperature transition from the disproportionated to a 
homogeneous phase is accompanied by vanishing of the charge gap. With the increasing crystal-field 
splitting the temperature range of the disproportionated phase shrinks and eventually disappears 
completely. 

PACS numbers: 71.10.Fd,75.30.Wx,71.30.-|-h,71.28.-|-d 
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The pressure or thermally driven spin state transitions 
play an important role in the physics of magnetic ox- 
ides [J. A notorious example is LaCoOa. Its peculiar 
magnetic and transport properties have attracted atten- 
tion of physicists for decades, yet interpretation of its 
behavior remains controversial. The main characteristics 
of LaCoOa are the temperature (T) dependencies of the 
magnetic susceptibility and the conductivity 0, Q , which 
exhibit three distinct regions: (i) a low-T non-magnetic 
insulator, (ii) an intermediate-T paramagnetic insulator 
and (iii) a high-T paramagnetic bad metal. It is gener- 
ally believed that the evolution from the non-magnetic 
to the paramagnetic state is due to thermal population 
of an excited state of the Co ion with spin (or spin and 
orbital) degeneracy. Commonly considered scenarios in- 
volve either low-spin to high-spin excitation 0, Q-Hi low- 
spin to intermediate-spin excitation or both @. 
Temperature affects the electronic system also indirectly, 
by changing the crystal-field splitting through the lat- 
tice thermal expansion. Yet another possible piece to 
the puzzle are deformations of CoOg octahedra and their 
coupling to the spin states of the Co ion [2], Q ■ Therefore 
it is rather difficult to distinguish the leading effects from 
the secondary ones. 

In this Letter we study a minimal fermionic lattice 
model that exhibits the spin state transition. As purely 
electronic it does not include the effect of lattice thermal 
expansion or the magneto-elastic coupling. The model 
is essentially the same as the one studied by Werner and 
Millis who used the dynamical mean- field theory 11 1 
(DMFT) to map out its phase diagram at fixed temper- 
ature, and by Suzuki et al. [l2j . who studied its ground 
state as a function of doping by variational Monte-Carlo. 
We employ the DMFT method to study the temperature 
dependent properties in the vicinity of the boundary be- 
tween the low-spin and high-spin phases. Unlike Ref. [13 
we assume a specific lattice, which allows us to investi- 
gate the ordering tendencies. We compute [l^ the one- 



particle spectra and the local as well as the uniform spin 
susceptibility and find that, similar to the behavior of 
LaCoOa, the model exhibits three distinct temperature 
regions. In addition to the low-T non-magnetic insulator 
and high-T local moment metal we find a spin dispropor- 
tionated insulating phase at intermediate temperatures. 
In order to interpret the DMFT results we construct an 
effective low-energy Potts model, which allows analytic 
calculations. 

Our starting point is a Hubbard Hamiltonian on a 
square lattice 
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where ^, h\ ^ ifli^m ^i.o-) are the fermionic creation (an- 
nihilation) operators for a spin index a and two types of 
orbitals for each lattice site i, and n^^, n\ ^ are the cor- 
responding occupation number operators. The nearest- 
neighbor hoppings taa = 0.45 eV, tbb ~ 0.05 eV and the 
on-site interaction parameters U = A eV and J = 1 eV 
are chosen to yield a broad a-band and a narrow b- 
band, mimicking the electronic structure of LaCoOa. 
The crystal-field splitting is denoted with A. The T- 
dependent chemical potential is fixed to yield the av- 
erage filling of 2 electrons per lattice site. Unlike Ref. [l3 
we consider only Ising terms in the on-site interaction. 

The context of the present study is set by a conceptual 
form [1^ of the [/ — A phase diagram of Ref. [13 shown 
in Fig. [TJ The boundary of the metallic phase is given by 
opening of a linearly increasing charge gap (indicated by 
the color intensity). The line separating high-spin (HS) 
Mott and low-spin (LS) band insulator corresponds to 
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FIG. 1: (Color online) Conceptual phase diagram of the two 
band model for (7/ J = 4. The size of the charge gap is indi- 
cated by color intensity (white=no gap). The parameters of 
the present study are marked with the cross. 



degeneracy of local HS and LS states. The parameter 
range of interest corresponding to small gap LS insulator 
is close to the triple point. In the following we present 
the results for A of 3.40 and 3.42 cV (marked with the 
cross in Fig. [T]). 

The DMFT equations are solved as described in Ref.fisl 
using the strong coupling continuous time quantum 
Monte-Carlo solver [16|. Guided by proposals of spin 
disproportionation in LaCoOs 0, [13 we have doubled 
the unit cell to allow for a spontaneous two-sublattice 
order (the sublattices are denoted with A and B). In the 
lower panel of Fig. [5] we show the occupation numbers 
n\ gas a function of temperature. Below 500 K a dis- 
proportionation takes place. The local states on the A 
sites can be described as HS-I-LS statistical mixture of 
LS (a°6^) and HS (a^b^) configurations with short excur- 
sions to 1-electron (a°6^) configuration. The B sites host 
the LS (a" 6^) states which experience short excursions to 
3-electron (a^b^) configurations. The driving force of the 
disproportionation is a superexchange interaction medi- 
ated by these excursions, as shown in the inset of the 
lower panel of Fig. [51 

In the upper panel of Fig. [2] we show the local sus- 
ceptibility averaged over A and B sites (the A and B 
contributions are shown in the inset). The local suscep- 
tibilities reflect the probability to find a given site in the 
HS state, while the occupancies fi^ g are to some extent 
affected by the charge fluctuations to 1- and 3-electrons 
states on A and B sites respectively. Comparison to the 
local susceptibility calculated in the (by constraint) ho- 
mogeneous phase reveals an enhancement of the average 
HS abundance due to the disproportionation. We have 
calculated also the uniform susceptibility by adding a 
small Zeeman field to Eq. [TJ In the disproportionated 
phase the uniform susceptibility coincides with the av- 
eraged local susceptibility. This is a simple consequence 
of the local moments on A sites being separated from 
each other by the B sites hosting the LS singlets. In the 
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FIG. 2: (Color onhne) The upper panel shows the T depen- 
dence of the average local susceptibility per atom for the crys- 
tal field parameter A of 3.40 (red) and 3.42 (black) eV. (The 
site resolved contributions are shown in the inset.) The data 
for A = 3.40 eV are compared to the uniform susceptibility 
(diamonds) and the local susceptibility obtained in the ho- 
mogeneous phase (empty circles). The lower panel shows the 
difference n'X — n% of the occupancies between the A and B 
sites. The dotted lines are guides to the eye. 



high-T homogeneous phase the uniform susceptibility is 
found to be enhanced over its local counterpart. This is 
somewhat counterintuitive since a naive expectation of 
an anti-ferromagnetic superexchange between the neigh- 
boring HS excitations should lead to an opposite effect. 
The Pauli susceptibility associated with the bad metallic 
state, discussed next, is an order of magnitude too small 
to provide an explanation. Therefore we conclude that 
an effective ferromagnetic coupling exists in the high-T 
phase. 

The evolution of the one-particle spectra is shown in 
Fig. [31 The disproportionated phase exhibits a well- 
defined charge gap which starts to fill with incoherent 
excitations as the system approaches the transition to 
the homogeneous phase. To interpret this behavior we 
consider the definition of the charge gap in terms of the 
eigencnergies of the system Ejsi+i -\- Ej^-i — 2En, where 
En corresponds to an eigenstate with a non-vanishing 
occupancy and i?7v+i, E^-i arc energies of lowest states 
that can be reached by adding or removing an elec- 
tron. At zero temperature the ground state is a prod- 
uct of LS [aPb"^) configurations on each site. The lowest 
iV + 1 state corresponds to a single a-electron propagat- 
ing over the the filled b band, while the lowest — 1 
state corresponds to a single 6-hole. The gap is obtained 
as the on-site contribution reduced by half-bandwidths 
t/ - 5 J + A - Wa/2 - Wb/2. 

The situation at elevated temperatures is more compli- 
cated as the initial states involve also sites in HS config- 
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FIG. 3: (Color online) The T evolution the one-particle spec- 
tra for A — 3.40. The spectral densities of different orbitals 
are resolved by color: a (black) and b (red), while the dashed 
lines correspond to the B sites and full lines to the A sites. 



uration as well as sites with one or three electrons. The 
disproportionatcd and the homogeneous phases differ in 
the constraints posed on the relaxation of -I- 1 and iV — 1 
excitations by the states of the neighboring sites. While 
in the homogeneous phase each site has neighbors which 
fluctuate into 1- and 3-electron states, in the dispropor- 
tionatcd phase fluctuations are either to 1-electron states 
on A sites or 3-electron states on B sites. As an example 
we discuss the fo-elcctron excitations on an A site from 
the HS initial conflguration. (The LS initial configura- 
tions contribute with completely filled valence band.) In 
the disproportionatcd phase the energy of lowest N + 1 
state reached by adding b electron consists of the on-site 
contribution reduced by Wa/2 due to a electron hopping. 
The lowest — 1 state reached by b electron removal cor- 
responds to (a^b^) conflguration on the A site and (a°6^) 
configuration on the neighboring B site. The correspond- 
ing gap estimate is still finite {U — 2J — Wa/2). In the ho- 
mogeneous phase the lowest iV — 1 excitation connected 
with b electron removal is the same as just mentioned. 
For -|- 1 excitations there is an additional possibility in 
the homogeneous phase to add b electron on a site in HS 
state while its neighbor is in (a"6^) configuration leading 
to the (a°&^) final state with {a^b^) on the neighboring 
site after a electron transfer. Considering these excita- 
tions we obtain a vanishing estimate for the charge gap. 

In order to gain insight into the DMFT results we in- 
tegrate out charge fiuctuations in ([T|) to get an effective 
Potts model with three lowest-energy states (LS, HSt, 
HSJ,) per site 

h-coy: + E i^^^^'^f' + ■ (2) 

Here nf^ and n^^ are the projectors of the three states. 



^2 = ~"U+j ■ should be pointed out that this model is 
only good for qualitative comparison to the DMFT data 
as the charge fiuctuations are not negligible (in particular 
in the metallic phase). A mean-field decoupling of ([2]) 
allowing for a two sublattice order leads to the free energy 
per site expression 

= y (-^A + Xb) + 2^1 (XA + XB - 2xaXb) - ilXAXB 

T T 
+ -(1 - xa) ln(l - xa) + -(1 - a^s) ln(l - xb) 



T XA. T XB. 

-XA In(^) + -^XB M-^)> 



(3) 



where xa,b are the mean values of n^^ + n^^ on the two 
sublattices. The equilibrium values of xa — xb obtained 
by minimization of ([3]) are shown in Fig. |4] together with 
corresponding uniform spin susceptibility. For > 
we find a uniform LS ground state at T — 0, which is 
followed by a transition into the disproportionatcd phase 
characterized by xa xb between temperatures Td and 
Tci- With increasing the Tc\ and Tc2 converge and the 
disproportionatcd phase eventually disappears for large 
enough ^o- 

The disproportionatcd phase exhibits an enhanced sus- 
ceptibility, which has two sources. First, like in the 
DMFT results the HS abundance is enhanced with re- 
spect to the homogeneous phase. Second, in the ho- 
mogeneous phase the anti-ferromagnetic (AFM) coupling 
^2 reduces the uniform susceptibility.Note that this is in 
contrast to the DMFT results. 

The HS-LS model is not new. It was suggested for 
LaCoOa by Goodenough and Raccah Q and the dis- 
proportionation was studied by Bari and Sivardiere [l^ ■ 
However, traditionally the interaction between the LS 
and HS state was associated with the magneto-elastic 
coupling, namely a breathing distortion of the CoOe oc- 
tahedra. This is quite different from the present work 
where the disproportionation is of purely electronic ori- 
gin. 

Finally, we discuss the implications of our results in the 
context of LaCoOa. We show that physics of LaCoOs 
over the whole temperature range can be qualitatively 
described by purely electronic effects and with only one 
local moment state. The simplicity of the present model 
lends it generality, but inevitably involves approxima- 
tions. First, our model involves two non-degenerate or- 
bitals while LaCoOs is characterized by three-fold quasi- 
degenerate t2g and two-fold degenerate Cg orbitals lead- 
ing to different entropic contributions. Nevertheless, the 
main control parameters are the coupling constants, and 
different band degeneracies have only quantitative effect. 
Second, the Co-0 bond-lengths change due to the nor- 



4 




FIG. 4: (Color online) The uniform susceptibility of the Potts 
model as a function of temperature and the parameter ^o- 
(upper panel) The difference of the HS populations xa — xb 
as a function of temperature and ^o- (lower panel) 

mal and anomalous lattice thermal expansion Q. The 
latter is caused by populating the HS state, well-known 
to weaken the metal ion- ligand bonds [l9|, [2^. Con- 
sidering a path of decreasing ^o{T) in Fig. 2] it is clear 
that the lattice response enhances the observed effects 
as smaller crystal field favors the disproportionation as 
well as the metallization. Third, no disproportionation 
and corresponding breathing lattice distortion was ob- 
served in recent experiments [^,[2l|. It is well known that 
mean-field approximations overestimate ordering tenden- 
cies and that a long-range order in the mean- field solution 
is often indicative of short-range correlations in the sys- 
tem. Therefore we speculate that in LaCoOa dynamical 
HS-LS correlations take place. It is plausible that such 
a dynamical effect can arise from the instantaneous HS- 
LS interaction of electronic origin. On the other hand, 
magneto-elastic HS-LS interaction retardation effects due 
to the lattice dynamics arc likely to weaken dynamical 
HS-LS correlation considerably. Fourth, two regions of 
Curie- Weiss behavior observed in experiments have been 
interpreted as an evidence for two different local moment 
states (high-spin and intermediate-spin) . In our model 
we observe different behaviors of the uniform suscepti- 
bility in the disproportionated and in the homogeneous 
high T phase. The difference has three sources: (i) an 
enhanced abundance of the HS configurations in the dis- 
proportionated phase (over a hypothetical homogeneous 
phase at the same T), (ii) absence of the nearest- neighbor 



anti-ferromagnetic correlations in the disproportionated 
phase, (iii) the metallicity of the homogeneous phase. 

In conclusion, we have used numerical DMFT method 
to study a two-band Hubbard model with quasi- 
degenerate high-spin and low-spin local states. Varying 
temperature we have observed three different regimes: a 
low-spin insulator, an insulating phase with HS-LS dis- 
proportionation and enhanced Curie- Weiss susceptibility 
and a homogeneous metallic phase with Curie- Weiss sus- 
ceptibility. We have argued that our model study cap- 
tures the essential physics of LaCoOa and thus that the 
properties of LaCoOs can be explained with a single mag- 
netic moment carrying state and without the effect of the 
lattice thermal expansion. 
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